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Abstract. We consider the existence and asymptotics for the fundamental solution of an 
elliptic operator in nondivergence form, C(x,d x ) = a,ij(x)didi, for n > 3. We assume that the 
coefficients have modulus of continuity satisfying the square Dini condition. For fixed y, we 
construct a solution of CZ y (x) = for < \x — y\ < e with explicit leading order term which is 
0(\x — y\ 2 ~ n e I ^ x ' v ^) as x — > y, where I(x, y) is given by an integral and plays an important role 
for the fundamental solution: if I(x, y) approaches a finite limit as x — > y, then we can solve 
C(x, d x )F(x, y) = 8(x — y), and F(x, y) is asymptotic as x — » y to the fundamental solution for 
the constant coefficient operator C(y,d x )- On the other hand, if I(x,y) — > — oo as x — > y then 
the solution Z y (x) violates the "extended maximum principle" of Gilbarg & Serrin [5] and is 
a distributional solution of C(x, d x )Z y (x) = for \x — y\ < e although Z y is not even bounded 
as x — > y. 

1. Introduction 

Background. Consider an elliptic operator in nondivergence form 

(1) £(x, d x ) u(x) = aij(x) did.ju(x), 

where di = d/dxi and we have used the summation convention for repeated indices. The coeffi- 
cients aij — ciji are real- valued functions defined on R" for n > 3, and we denote the symmetric 
and positive definite matrix (cnj(x)) by A x . (The case n = 2 can be treated with a similar 
analysis, but additional complications arise which we have chosen to avoid here.) A fundamental 
solution for £ in an open set U is a function F(x,y) satisfying F(x, •) £ Lg oc (U) and 

(2) -£(x 1 d x )F(x 1 y) = 5(x-y) for x, yeU 

in a distributional sense that needs to be made clear; for this some regularity of the coefficients 
will be required. If F(x, y) satisfying ^ exists, then the operator Q defined by 



(3) Q<P(x) = ~ / F(x, y) 4>(y) dy 

Ju 

provides a right-inverse for £ on Co(U), the space of continuous functions with compact support 
in U. 

In the "classical" case that the coefficient functions are A-H61der continuous in a bounded 
domain U for some A S (0, 1), it is well-known (cf.[21]) that a fundamental solution exists in U 
and is asymptotic (as x — > y) to the fundamental solution for the constant coefficient operator 
obtained by freezing the coefficients ay at y: for n > 3 this means 

(4) F(x,y) = F y (x-y)(l+H(x,y)), 
where, letting (,) denote the inner product in R", 

rx\ p i ^ (Ay 1 x, x) ^ 

(5) FJx) = 

W yy ' (n-2)\S n - 1 \^/MA^ 
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is the fundamental solution for the constant coefficient operator C(y,d x ) — aij{y)didj, and the 
remainder term H(x,y) in ([4]) satisfies 

(6) \H(x,y)\+r\D x H(x,y)\+r 2 \D 2 x H(x,y)\<cr x as r = \x - y\ -» 0, 

for all y in a compact subset of U. This fundamental solution may be used (cf.[3T]) to show 
the classical regularity result: if u € C 2 (U) and Cu is A-H61der continuous in U, then didjU is 
A-H61der continuous in U. 

The Holder continuity may be generalized by assuming the coefficients have a weaker modulus 
of continuity, i.e. a»j € C U (U) where w(r) is a continuous, nondecreasing function for < r < 1 
satisfying oj(Q) = 0, and 

(7) C"(U) = {feC(U):\f(x)-f(y)\<coj(\x-y\) for x, yet/}. 

If lu satisfies the Dini condition at zero, i.e. w(<)t _1 dt < oo, then we say that the coefficients 
are Dini continuous. In this case, there are regularity results analogous to the case of Holder 
continuity (cf. Proposition 1.14 in Chapter 3 of [23]); however, we could not find in the literature 
an asymptotic description of the fundamental solution such as Q with estimates on the second- 
order derivatives D 2 H(x,y). 

Dini continuity is also essential for the "extended maximum principle" of Gilbarg & Serrin [8] 
to hold: a C 2 -solution of 

Cu > for < \x\ < r 

with 

u(x) = o{\x\ 2 - n ) &s\x\ ^0 

must satisfy 

u(x) < M := max{u(y) : \y\ = r^} for < \x\ < ro, 

and limsupi^i^Q u(x) < M. In fact, they give an example (which we will discuss in Section 2) in 
which the coefficients are not Dini continuous and the extended maximum principle fails. 

The above regularity assumptions (Holder or Dini continuity) on the coefficients are required 
to study the behavior of the fundamental solution as a function of x (for fixed y) . If we instead fix 
x and consider the behavior in y, then regularity of the coefficients is not required; however, 
we cannot expect to achieve as precise an asymptotic description as (JU). This is most conveniently 
described in terms of the Green's function for ([1]) on a smooth, bounded domain U, which may 
be defined (as in [3]) to be G(x, •) e L\ oc (U) satisfying 

(8) 4>{x) = - / G(x, y) C(y, d y ) <j>(y) dy for S C 2 (U) with = on dU. 

Notice that ([5]) can be expressed formally as —£* (y, d y )G(x,y) — 5(x — y) and implies that 
Q4>(x) = — Jjj G(x, y) 4>(y) dy defines a left-inverse for C on C 2 (U). When the ay are measurable, 
bounded, and uniformly elliptic in U, then the Green's function is known to exist and Fabes 
& Strook [6] showed that G(x, ■) S L q (U) for some q > n/(n — 1), while Bauman [I], [2], [3], 
and Escauriaza 4J obtained pointwise estimates on G(x,y) as y — > x in terms of a nonnegative 
"adjoint solution" W(y) which satisfies £* (y , d y )W (y) = in U. However, our paper is not 
concerned with such general coefficients, and for us a Green's function will also be a fundamental 
solution in the sense of 

Our results. In this paper, we allow the coefficients to be less regular than Dini continuous, 
and we want to study the solutions of Cu(x) = with an isolated singularity at x = y, as well as 
the existence and asymptotics of a fundamental solution F(x, y) satisfying ^ in an appropriate 
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distributional sense. We assume that the coefficients have modulus of continuity ui satisfying the 
"square-Dini condition" 



(9) 



f 1 2, x dt 

/ LO Z {t) — < 00. 

Jo t 



Condition ((9]) has been encountered by other authors in different contexts: cf. [5], [7], [IT], [23] . 
To construct our solution of ((2]), we first fix y and seek a solution of 

(10) C{x,d x )Z y (x)=0 for x e B e (y)\{y}, 

where B £ {y) — {x : \x — y\ < e} for e sufficiently small, and Z y (x) has the appropriate singularity 
as x — » y. Assuming that the modulus of continuity at y satisfies (j9)), we shall construct a solution 
of ([TO]) with the asymptotic description 

(11) Z y {x)~ (A- 1 (x-y),(x-y)) 2 -^e 1 ^^ as x -» j,, 

where the factor e 1 ^^' adjusts for lack of regularity in the coefficients: if the <Xij are Holder 
continuous, then we can take I(x,y) = and c y Z y (x) is asymptotic to F y (x — y) as x — > y. In 
general, however, we find that 



(12) /(x,y)=/ J/ ^ v /(A I 7 1 (a ; -y),( : r-y)) 
where I y (r) is given by 

(13) I ( tr^A" 1 ) - n .< W /S (* - A,7 1/2 (z - y)) ^ (b 



with tr denoting matrix trace. As r — > 0, / y (r) need not even be bounded, so the singularity of 
Z y (x) need not be 0(|x — yj 2- ") as it was in the Holder case. These formulas simplify significantly 
if we use an affine change of variables in which y corresponds to x — and (0) = S^j : 

(14) Z{x) - \x\ 2 - n e I{lxl) as \x\ -> 0, 
where 

(15) 7(r) = Is^i / <N< , ( tr<Al) " " S?/ FF' 

We can verify that the absolute value of the integrand in (fT5|) is bounded by w(|z|), so the coeffi- 
cients being Dini continuous implies that the improper integral defining 1(0) converges absolutely. 
Even if 1(0) does not converge, we shall see that for any A > there exists C\ > such that 

(16) \I(r)\ < A|logr| +C X for < r < e, 

so the singularity of Z at x = is never very far from |a:| 2_rl . Nevertheless, the behavior of I(r) 
as r — > plays an important role for our results. There are three principal cases to consider. 

1. 1(0) = lmv—^o I (t) exists and is finite. 

In this case, we can scale Z(x) by a constant multiple to obtain a solution that is asymptotic to 
the fundamental solution for the Laplacian. In fact, we can solve the distributional equation 

(17) - C(x,d x )Z(x) = C 6(x), 
and find 

(18) C Q = (n-2)\S n - 1 \e I{0 \ 

Note that the improper integral defining 1(0) may converge even if the modulus of continuity 
does not satisfy the Dini condition. 
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2. I(r) — > — oo as r — > 0. 

We see that Z(x) = o(|a;| 2 ~ n ) as |x| — ► 0; we still have Z[x) — ► +oo as x — > 0, so this violates the 
extended maximum principle of [5] . Nevertheless, we can solve (fTT|) to find Co = and this yields 
the surprising fact that there exists a distributional solution of C(x,d x )Z(x) = for x £ B E {y) 
which is not even bounded. 

3. I(r) — > oo as r — > 0. 

Now we find (x)|iz;|™ _2 — > oo as \x\ — » 0, so this solution grows more rapidly than the funda- 
mental solution for the Laplacian. However, we cannot solve (fl7|) for Co in this case. 

We next allow y to vary over U. Provided that we are in Case 1 at each y S U, we can use 
Z y (x) to construct the fundamental solution F(x,y) in t/. Our main result (Theorem 3) states 
that, provided € C"(U) where to satisfies © and I y (0) = lim r ^o ly i r ) exists and is finite 
for each y G U, a fundamental solution F(x,y) exists in the form ([4]) where the remainder term 
H(x,y) may be estimated in LP for any p S (l,oo) in terms of a? and the rate of convergence 
I y (r) — > /y(0). More specifically, let us assume u(r)r~ 1+K is nonincreasing for < r < 1 where 
k € (0, 1), and introduce 



where 6{r) is a nondecreasing function of r satisfying 6(0) = 0. Then for y in a compact subset 
of U, we have 



where M p (f(-,y),r;y) denotes the L p -mean of f(x,y) as a function of x (for fixed y) over the 
annulus A r (y) = {x : r < \x — y\ < 2r}. (In fact, when lu satisfies the Dini condition, the 
right hand side of (p?Tj) reduces to just cui(r).) Taking p > n, we obtain pointwise bounds on 
\H(x, y) \ + r\D x H(x, y)\, but we no longer have pointwise bounds on D^.H(x, y) as we did in © 
when the coefficients ay were Holder continuous. 

Note. In (|2Tj) and throughout this paper, c is used to denote a constant whose value may change 
line-by-line. It may depend upon n and the dy, but not on r. 

Organization of this paper and comparison with other works. The organization of this 
paper is as follows. In Section 1 we discuss some preliminary estimates for solutions of the Poisson 
equation Au = f. In Section 2 we construct the solution Z y (x) of (fiT)|) by first considering the 
case (0) = Sij and then performing a change of coordinates. In Section 3 we analyze the 
equation —C(x,d x )Z y (x) — C y S(x — y) and calculate C y when I y (0) is finite or — oo. Finally, in 
Section 4 we apply the previous results to construct our fundamental solution F(x, y) in the form 
([4]) and obtain estimates on the remainder term H(x,y). The analysis in each section makes use 
of Lp-means on annuli A r = {x : r < \x\ < 2r} to measure the growth and decay of functions as 
r — * 0; these are defined in Section 1. 

The results obtained and the techniques used in this paper are closely related to those in our 
previous paper [T7] . In [T7] we studied the asymptotics for solutions of the adjoint equation 
for C, and discussed their relationship to the estimates obtained in 0] as y — > x. In fact, the 
condition that I y (r) is bounded (above and below) as r — > not only allows one to conclude 
that the solution Z y of (jTUJ) satisfies c y \x — y\ 2 ~ n < Z y (x) < c' y \x — y\ 2 ~ n as x — > y, but 
that C* (x, d x )u(x) — admits solutions in B e (y) whose L p -mean is bounded between positive 
constants as r — \x ~ y\ — * 0; if I y (0) exists and is finite for every y, then this bounded solution 
of C* (x, d x )u(x) = is continuous. 



(19) 





\I y (r) - I y (0)\ < 6{r) for all ye U, 



(21) 



r 2 M P (D H (• , y) , r;y) < c max(w(r), cr(r), 6(r)) as r — > 0, 
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Both [17j and the present work are independent of, but related to, the asymptotic theory 
developed in [14]. In particular, L p -means were extensively used in [12] and [13] . The asymp- 
totic formulas that we obtain are analogous to those of [15] , where an asymptotic representation 
near the boundary was obtained for solutions to the Dirichlet problem for elliptic equations in 
divergence form with discontinuous coefficients. In particular, note that we do not make use of 
the maximum principle in this paper; in fact, most results of this paper hold for complex-valued 
coefficients, although this introduces some technical difficulties which we have chosen to avoid. 

2. Preliminaries 

Throughout this paper, to measure the growth of functions as x — > y, it will be convenient to 
use the L p -mean for some p £ (1, oo): 

(22) M p {w,r,y)=(p \w(x)f dx] , 

V Mv) j 

where A r (y) — {x : r < \x — y\ < 2r}; here (and elsewhere in this paper) the slashed integral 
denotes mean value. It is natural to also introduce 

(23) Moo(w,r;y) = sup \w(x)\. 

Mv) 

We may apply (I22| and ([23|) to vector or matrix valued functions w with \w\ denoting the norm. 
We also need to estimate derivatives, so we define 

(24) M 1<00 (w,r;y) = rM^Dw^y) +M ao {w,r;y), 
where Dw represents the gradient of w, and for p€ (1, oo) 

(25) M 2 , p (w, r- y) = r 2 M p {D 2 w, r; y) + rM p (Dw, r; y) + M p (w, r; y), 
where D 2 w represents the Hessian matrix of w. Sobolev estimates show that 

(26) p > n =>■ M lt00 (w,r;y) < cM 2 , p {w,r;y). 

When y = 0, we shall abbreviate M p (w,r; 0) as M p (w,r) (and similarly for Mi :OQ and M2, p ). 

For x € K™\{0}, let 9 = x/\x\ £ S'™~ 1 and let d9 denote the standard surface measure on 
S n ~ l , We will use the spherical mean of a function w. 



(27) W(r) = r w{r6) dO. 

In particular, in this section we consider the equation 

(28) Av = f in R"\{0} 

when / = 0, and investigate the behavior of the L p -mean of the solution as x — * 0; our results are 
quite analogous to those of [12] and [14]. We shall let F(|a;|) = c„ |x| 2_ ™ denote the fundamental 
solution for the Laplacian in R n . 

Proposition 1. Suppose n > 2, p S (l,oo), and f £ L^ oc (W\{0}) satisfies / = 0, 

(29) / |aj||/(x)|da;<oo ) and f lli^A. dx < oo. 



J\x\<l J\x\>l \ x \ 

Then v = K f = T * / defines a distribution solution of (28\) that satisfies 

(30) M 2 , p (Kf,r)<c(r 2 M p (f,r)+r 1 - n [ \x\ \f(x)\ dx + r / l /j^dx 

\ J\x\<r J\x\>r \ x \ , 
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where we have introduced 

M p (w,r) := (f \w(x)\ p dx] 

\J r/2<\x\<4r J 



Elementary estimates can be used to show 



(31) 



|a;|<r 



\g(x)\dx<c / M p (g,p)p n ~ 1 dp 



with an analogous estimate for \x\ > r. These estimates enable us to obtain the following corollary 
from Proposition [TJ 

Corollary 1. Under the conditions of Proposition^ 

(32) M %p (Kf, r) < c r 1 ^ / M p (f, p)p n dp + r M p (f, p) dp. 



Proof of Proposition 1. Let v = Kf and let \ denote the characteristic function for the 
annulus A r = B^ r \B r / 2 - Using / = 0, we see that 



v(x)= T(\x - y\)x(y)f(y)dy ■ 



(T(\x-y\)-T(\x\))f(y)dy 



\y\<r/2 



{T(\x-y\)-Y{\y\))f{y)dy. 



\y\>4r 



We want to estimate M p (v,r), rM p (Dv,r), and r 2 M p (D 2 v,r), and show that each is bounded 
by the right hand side of (|30|) . 

Using Stein's inequality [52], for < a < n/p and < b < n/p' with a + b = 2 we have 



^{\x-y\)x(y)f(y)dy 



< cr°||x(l/)|y| 6 /(i,)llzw 



LP(A r ) 
\b JV„,\II „„2| 



= cr a \\\y\ b f(y)\\ LP{Ar) <cr 2 \\f\\ LP{Ary 
It is elementary to show that for \y\ < r/2 and r < \x\ < 2r we have |r(|a; — y\) — r(|x|)| < 



c a; 



l-n 



y\, where c is independent of r, so for x € A r we have 



(T(\x-y\)-T(\x\))f(y)dy 



\y\<r/2 



< c\x 



l-n 



\y\<r/2 



\f(y)\\y\dy. 



Similarly, we can show that for \y\ > 4r and r < \x\ < 2r we have |r(|a; — — r(|y|)| < cl^llyl 1 ™, 
so x G A r implies 



(T(\x-y\)-T(\y\))f(y)dy 



<c\x\[ \lM-dy. 

J\y\>4r \y\ n 



Using these estimates, we easily conclude that M p (v, r) is bounded by the right hand side of ([ST))) . 
Next we consider 



d t v(x) = I Ti(x-y)x(y)f(y)dy+ / T l {x-y)f{y)dy+ f Ti(x - y)f{y) dy, 

J B r / 2 



where 



r i (x) = r'(\x\)-^- 
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Applying Stein's inequality as above but with a + b = 1, we conclude 



Fi(x-y)x(y)f(y) d v 

Elementary estimates for r < \x\ < 2r show that 



< cr I 



LP(A r ) 



^i{x-y)f{y)dy 



r/2 



(Ti(x-y) - T l {x))f{y)dy 



I LP(A r ) ' 



< c|x|" 



\f(y)\\y\dy 



\y\<r/2 



and 



Fi(x - y)f(y)dy 



(T i (x-y)-T i (y))f(y)dy 



< c 



U, in— 1 a 
|j/|>4r 12/1 



From these estimates we easily conclude that rM p (Dv,r) is bounded by the right hand side of 
(US). 

Finally, we consider 

didjv(x) = / Tij (x - y)x{y)f{y) dy + T l3 [x - y)f(y) dy + T l3 (x - y)f(y) dy, 



where is the singular kernel given by 



r^x) = r"(\x\& + r'(\x\)-^ x] 

\X\ 



Using the L p -boundedness of singular integral operators, we conclude 



Fij(x - y)x(y)f(y)dy 



<c\\f\\ 



LP(A r ) 



LP(A r )- 



Elementary estimates for r < \x\ < 2r show that 



Fij(x-y)f(y)dy 



< c\x\ 



{Tij(x -y)- Tij(x))f(y) dy 



\f(y)\\y\dy 



\y\<r/2 



and 



/ V i3 ix-y)f{y) 



dy 



{T t3 (x-y)-T l0 {\y\))f{y)dy 



<c\x\~ l j ^dy. 

J\y\>4r \y\ 



These estimates show that r 2 M p (D 2 v, r) is bounded by the right hand side of (|30[) . This completes 
the proof. □ 



3. Constructing the Singular Solution Z y (x) in B £ (y)\{y} 

In this section, we fix y and construct a solution of (flu]) for e sufficiently small. Since the result 
is localized near y, the ellipticity and continuity of the coefficients of L need only be assumed at 
y. In fact, we first consider the case when y = and a 2 j(0) = 5^: 



(33) 



sup || A x - I|| < uj{r) for < r < 1, 

\x\ — r 
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where I is the identity matrix and w is a continuous, nondecreasing function for < r < 1 
satisfying the square-Dini condition |9]). For convenience, we shall also assume the monotonicity 
condition that for some k G (0, 1) we have 

(34) oj(r) r~ 1+K is nonincreasing for < r < 1. 

The significance of Q34j) is that it requires w(r) to vanish more slowly than r as r — ► 0, which is a 
natural constraint. (As long as ui(r) r~ 1+K is nonincreasing for < r < e with some e > 0, then 
u> may be modified for e < r < 1 to satisfy (|3"4")1 .) We seek a solution of 

(35) C(x, d x ) Z(x) = for x £ B £ \{0}, 

where B e = {x : \x\ < e} with e sufficiently small. Although we generally assume that C has 
real-valued coefficients, the theorem below holds when the are complex- valued. 

Theorem 1. For n > 3 and p £ (l,oo), suppose the coefficients dij(x) are bounded, measurable 
functions satisfying \33\) . For e > sufficiently small, there exists a solution of 135\) in the form 

(36) Z(x) = h(\x\)+v(x), 
where h is of the form 

(37) h(r)= /"V- n e J «<i8(H-C(r)), 



wrai/i 7(r) given by 115\) and 

(38) M 2 , P (C, r) < c max(w(r), a{r)), 
where a is given in \19\) . and v in R36]) satisfies 

(39) M 2iP (v,r) < cr 2 - n e'( r > W (r). 

Moreover, for any u e W^(B £ \{0}) that is a strong solution of C{x,d x )u = in _B e \{0} subject 
to the growth condition 

(40) M 2 , p (u,r) < cr 1 "™ 4 " 60 where e > 0, 
i/iere eiisi constants C, Co, Ci, . . . , C„ (depending on u) such that 

n 

(41) u(af) 

3=1 

where w satisfies 

(42) M 2 ,p(w,r) < cr 2 ~ £l /or any e x > 0. 



We shall prove this theorem below, but first let us make some observations. In general, we 
do not know whether I(r) is bounded as r — > 0, but we can verify that < cr~ l Lo(r), so 

integration by parts in (|3"Tj) shows that 

r 2— n 

(43) /i(r) = -e^W + ^fr), 

n — z 

where /ii(r) satisfies Mi j00 (/ii,r) < cr 2_ "e J W max(u(r), cr(r)). If we take p > n and apply (f2"6"| 
to w, we conclude that 

I |2-n I(\x\) 

(44) Z{x) = ^ (1 + ^)) asM^O. 

n — I 

where Mi (0O (£,r) < c max(u>(r), cr(r)). Obviously, we can multiply the Z of (|4*4")) by n — 2 to 
obtain the Z of (IB 
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Even when I(r) is not bounded as r — ► 0, we can derive useful bounds on Z(x) as |ar| — > 0. It 
is not difficult to verify that a symmetric matrix A satisfies 

(45) -2(n-l)||A-I|| <tr(A)-n(Ay,y)|yr 2 <2(n-l)||A-I|| for \y\ = 1, 
so there exist constants c, C > so that Z satisfies (assuming n > 3) 

(46) c |af "™ exp (^-c n w(t)*J < |Z(x)| < C |a;| s 

as a; — > 0, where c„ = 2(n — l)/|S f,l ~ 1 |. Using (|46p and the fact that w(r) — > as r — > 0, we obtain 
(|16p . which shows that the singularity of Z at x — is very close to the classical case. 
An interesting class of examples is obtained by letting 

(47) a tj (x) = 5ij + g(\x\)xiXj\x\~ 2 , 

where g(0) — but vanishes slowly as r — > 0. Gilbarg & Serrin ,8 a used (|47|) with certain specific 
functions g to show that Dini continuity is essential for their extended maximum principle to 
hold. In our formulation, 

(48) "(A,)-»^ = (l-n)#|), 



dp 



so 

(49) I(r) = (l-n) [ S g(p) 

Jr P 

Thus any g(r) > which does not satisfy the Dini condition at r = (but does satisfy the square- 
Dini condition) will yield I(r) — > — oo as r — > 0, so the Z(x) of Theorem 1 gives an example of 
a solution of ([35]) with singularity at x — even though Z(x) = o(|x| 2- ™) as |x| — * 0, i.e. the 
extended maximum principle fails; the specific function in [8] is g(r) = —(1 + (n — 1) logr) . 

Proof of Theorem 1. Instead of showing the existence of Z in a very small ball B e , we shall 
replace the condition that u> satisfies the square-Dini condition by 



(50) 



where S is sufficiently small, and show existence in the unit ball B\. In fact, using (f50|) and (|34|) . 
we see that 

(51) w(r) < c K}Tl V5 for < r < 1, 
where c Kjn depends only on k and n: 

«> T ^{t) d ±>J{r)r~^ f t^dt = ^(r) c' K>n . 

Jr/2 1 Jr/2 

Moreover, it will be useful to consider C as defined on all of W L with C — A outside of B\. 
Therefore, we shall assume that 

(52) a ij( x ) — Sij for |x| > 1, 

and investigate a solution of CZ = in M™\{0}. 

To construct Z(x), we let h(r) = Z[r) denote the spherical mean as in |27|) . and let v(x) = 
Z(x) — hi\x\), so v(r) = 0. We shall reduce the problem to solving an operator equation of the 
form (I + S + T)v = f where S and T have small operator norm on a Banach space X defined as 
follows: for fixed p £ (1, oo), let us consider the functions v in 

W £oc for which the 

norm 

(53) \\v\\x= sup " — +sup £^-7= 

o<r<i ^(rje 1 ^' r >i v5 
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is finite, and take the closure to form a Banach space X. The reduction to the operator equation 
will take several steps. 

If we take the spherical mean of C(h + v) = 0, we obtain an ordinary differential equation for 

h: 



a {r) h » + a " (r) a(r V + aijdidMr) = 0, 



where 
(54) 



a n{ r ) — Jf a u{ r &) dO and a{r) 



S n-1 



Uij {r6)0i6jd0. 



From ([55)1 it is clear that these functions satisfy 

(55) l a ra( r ) — n \ < w ( r ) an d W{r) — 1| < w(r) for < r < 1. 

Using the fact that Av = Av — 0, the equation for h becomes 

(56) 



K + i,,(),i),rir: = 0, 



where @ij(x) = a,ij(x) — 5ij satisfies < ^d^D- Notice that v satisfies the partial differential 

equation 

(57) 



- Av = Pijdidjh - (3i j did jh + (iijdidjV - (3^8,8^. 

We shall first solve (|56|) and use that to eliminate h from (|57|) ; then we will be able to apply the 
results of Section 1. 

Let us introduce q = h! and rewrite (1561) as 



(58) 
where 
(59) 
and 



9 ' + n - 1 + R{r) 9 = B[D*v}(r), 



a n (r) 
R(r) = — — — n 



a(r) 



Notice that 

(60) 

and 



B[D 2 v](r) = -a~ 1 {r)/3 lJ d l d J v(r). 
\R(r)\ < coj(r) for < r < 1 and R(r) = for r > 1, 



\B[D 2 v](r)\ < CLu{r)\D 2 v(r)\ for < r < 1 and B[D 2 v](r) = for r > 1. 
Moreover, the monotonicity of ui{r) together with (|34[) imply 

(61) 

so we consequently obtain 

(62) M p (B[D 2 v} : r) <cuj(r)M p (D 2 v,r) for < r < 1. 

To solve (|55|) , let us introduce 



max Lo(p) < cuo{r), 

r<p<2r 



(63) 



E± (r) = exp 



r°° dt~\ r f 1 

± / R(t)— =exp ± / R(t) 

J r ^ _ . Jf 



dt 

T 



i 



E T (r) 



and observe that E±(r) = 1 for r > 1. It is useful to observe that 

dt 
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so as a consequence of (|60|) and (j5TJ) , we obtain 



(64) 



< p < r < 1. 



In particular, 

(65) ci E± (r) < E± (p) < c 2 £7± (r) for r < p < 2r, 
and for any g £ L^ oc (R"\{0}) we can readily verify that 

(66) M p (\x\ v E ± (\x\)g(x),r) < cr» E ± (r) M p (g,r), 
for any fixed 

It will be more convenient for us to use E±(r) than e ±/< - r - ) , but these functions are equivalent: 
if we note that (fTS"]) can be written as 

(67) I(r)= [\a n (p)-na(p)}^, 

Jr P 

then we see that 

(68) E+(r) = 4e ; W(l+r(r)), 
where A — exp^ 1 R(p)[l — a(p)]p~ 1 dp] is finite and positive, and 

dp 



(69) 



r(r) = exp 



R(p)(l-a(p)) 



satisfies |r(r)| < c<r(r). Thus for some constants c\,C2 we have 

(70) ciE + (r) < e I{r) < c 2 E + (r) for < r < 1. 
Now if we introduce </>(r) = r n ~ 1 E-(r)g(r), then we can rewrite (|58|) as 

(71) 0'(r) =r"- 1 £_(r) J B[ J D 2 u](r). 
But (fTTj) may be integrated to obtain 

(72) 0(r) = 0(0) + f p"- 1 ^ (p)B[£> 2 «] (p) dp, 



where 0(0) is an arbitrary constant. Of course, to conclude (|T2"|) , we must verify that 0' is 
integrable on (0, 1). But v £ X implies M p (D 2 v,r) < coj(r)r~ n E + (r), so we can use (|6"Tj) . (|6"2")l . 



(|65p . and Holder's inequality to obtain 
^ p n - l E_{p)B[D*v]{p)dp 

J r 



<cE_{r)u{r) / |D 2 w(a;)|dx < c w 2 (r). 

r<|x|<2r 



Now if we write 



/ p n - x E-{p)B[El l v]{p)dp = Y t p n - l E.(p)B[D 2 v](p)dp, 



then we obtain the estimate 



(73) 



/ p n -'E.(p)B[D 2 v](p)dp < c J2^(At) 



< c / (jj 2 (p) — = ccr(s) < cS. 
P 
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In particular, we see that / is integrable and so (|72|) is valid. This enables us to write 



h'(r) =g(r) = r^E+ir) 



(74) 
and 

(75) h"(r) = 1 - n - R ^ E +{ r) 



0(0) 



0(0) 



+ [ p n - x E_(p)B[D*v]{ P )dp 
Jo 

p n - 1 E^(p)B[D 2 v}(p)dp +B[D 2 ](r 



We can now use ([74]) and (|75|) to express 



f3 ij d i d j h = r- n E+(r) 



0(0)+ / p^E^B^v^dp 
Jo 

+B[D 2 v]{r)l3 t] 9 i e j , 



where 

(76) f{r6) = p u (r6) - (n + 11(^)13^6)6^ 

also satisfies |^(r0)| < ccj(r). Thus we can apply A" = T* to (f57|) to obtain 

(77) v + Sv + Tv = 0(O)w, 
where 

w(x) =K y ^ x [\y\- n E + (\y\)W(y)-^(\y\)}, 
\y\- n E+{\y\) 1* p^E^BiD^dp^iy)-^) 



Sv(x) = -K v ^ x 



and 



Tv = -K [B[D 2 v] (Pijdidj - fcOiO j) + fyjdidjV - fojdidjv] . 

For a given value of 0(0), we can solve ([77]) provided we show that w £ X and both 5 : X — > X 
and T : X — > X have small operator norms. 

To show w E X, we must estimate M^, p (w, r) for < r < 1 and for r > 1. For < r < 1 apply 
Proposition 2 to f(x) — \x\~ n E+(\x\)(ip(x) — ip(\x\)) (which vanishes for |x| > 1) to conclude 



M a , p (t«, r) < c \ r l ~ n J E+(p)uj(p) dp + rj E+(p)cj(p)p- n dp 
Using (l64|) and the fact that r 1 ~ K io(r) is nondecreasing (since both r 1 ~ K and uj are), we find 

E+(p)w(p) dp < E + (r) r c ^ 5 f P ~ c ^uj{p) dp 



< E+(r) io(r) r 



•-VS+1-K 



Taking S small enough that K — \ c5 > 0, we obtain 



p -cVs-i+K dp _ 



(78) 



E + (p)uj(p) dp < cr E + (r) u)(r). 



Using (|64|) and the fact that r 1+K uj(r) is nonincreasing (by ([34]) ). we find 



E + (p)Lo{p)p- n dp<E+{r)r- 



P C ^- n L0{p)dp 



< E + (r) ui(r) r 



p cVS+l- K -n d(J 
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For 5 small enough that n — 2 + k — c\fE > 0, we have 

(79) [ E + (pMp) p- n dp < cr 1 -" E + (r) w(r). 

J r 

Using |[78j) and |79}, we obtain 

M2 ^ r) :T < c forall 0<r<l. 
uj(r)E+{r) 

We can then use (|T0[) to replace E + (r) by e 1 ^ as required in the norm for X in Meanwhile, 
for r > 1 we use uj(1) < cVd and E + (p) < p~ cv ^ for < p < 1 to find 

M a , p (tu,r) < cr 1 -™ / E + {p)u(j>) dp < cv^r 1 "" p~ cVI dp < c \/£ r 1 "™, 



o 

provided J is sufficiently small. Consequently, 

— = < c for all r > 1, 

Vs 

and this confirms that w G X. 

Next let us show that S maps X to itself with small operator norm. We suppose that \\v\\x < 1 
and estimate M 2 p (Sv,r) separately for < r < 1 and for r > 1. For < r < f, we examine 
the proof of ([75)1 and observe that the condition \\v\\x < 1 enables us to choose the constant c 
independent of v. Thus the function 

fx(y) = \v\- n E+{\y\) ( V p^E-WBl&vMdpMy)-^])) 
Jo 

satisfies 

M p (/i,r) < c5E + {r)uj{r)r- n for < r < 1 
and M p (/i, r) = for r > 1 . For < r < I , we apply Proposition 2 to Sv — —Kfi to obtain 

M 2tP (Sv,r) < cS fr 1 -" f E + {p) uj{p) dp + r [ E + (p) u{p) p~ n dp) , 



and then use (J75J) and (|T9"|) to conclude (for S sufficiently small) 

M ^ r >^<cS forall 0<r<l. 

On the other hand, for r > 1, Proposition 2 implies (for (5 sufficiently small) 

M 2tP {Sv 1 r) ^cSr 1 - 11 E+{p) u(p) dp < c6 3/2 r 1 "" / p~ cv ^ dp < c <S 3/2 r 1 "". 
Jo Jo 

Thus we have 

M 2 , p (^,r)r"- 1 

— = < co for all r > 1. 

Combining these inequalities, we see that S : X — ► X has small operator norm. 

Finally, we show that T maps X to itself with small operator norm. We suppose that ||u||x < 1 
and estimate M2. p (Tv,r) separately for < r < 1 and for r > 1. Notice that the function 

h = B[D 2 v\ i t./jjK i7^r<) 

satisfies 

M p (f 2 ,r) < uj{r)M p (B[D 2 v],r) < cuj 3 {r) E+(r) r - n for < r < I , 
where c is independent of v, and M p {f 2 , r) = for r > 1 . Similarly, the function 
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satisfies 

M p (f 3> r) < cu{r)M p (D 2 v,r) < oj 2 {r) E+(r) r~ n for < r < 1, 

and M p (f 3 ,r) = for r > 1. For < r < 1, we apply Proposition 2 to Tv = -K{f 2 + fs) to 
obtain 



M 2 , p (T«,r) <c r 1 -" / a/(p) E+(p) dp + r oj\p) E + {p)p~ n dp . 



Using (5U, (El, and © 



M 2>p (Tu,r) r" 



< cv 7 ^ for all < r < 1. 



w(r)E+(r) 

On the other hand, for r > 1, we use (|5ip and l|64p to estimate 



M 2 , p {Tv,r) <cr l ~ n / £+(p) dp < c^r 1 



p- cV ~ s dp KcSr 1 



Consequently, 



M 2tP (Tv,r) r^ 1 



< c VS for all r > 1 . 



Combining these estimates, we see that T : X — » X has small operator norm. 

Since both £ and T have small operator norms on X, we conclude that (|77|) has a unique 
solution u, depending on the choice of the constant c\ = (f)(0). But once c\ and w are known, we 
obtain g(r) from (|74[) , and h(r) by integration of g(r): 



(80) 



/i(r) = 



ci+ / p^E-^B^v^dp 
Jo 



ds + C2 . 



where c 2 is an arbitrary constant. To obtain the desired solution of Theorem 1, we choose C\ 
to enable us to replace E + (r) by e 1 ^ for < r < 1. Using ([58"]) we see that we should choose 
c\ = A^ 1 and write h(r) = h (r) + hi(r) + c where 

(81) h (r) = [ s^V^ds 
and (recalling r from 

(82) hi(r)= [ s 1 ~ n e I ^T(s)ds+ f s l - n E + {s) ( p"' 1 E_(p)B[D 2 v}(p)dpds. 

Jr Jr JO 

Now integrate by parts to obtain 



ho(r) 



„2-n 



+ C + 



1 



n - 2 n - 2 

But |/'(s)| < cu>(s)/s and, similar to (|64|) . we can show that 



2 - n e /(s) /'(s)ds. 



zVS 



for s > r, 



so we may use to obtain 



< CCJ ( r ) e ^(^) r -^-l+«[ r 3-n+cV5-K + !] < cr 2-n w ( r ) e /(0 

provided <5 is sufficiently small. Thus we find 



(83) 



„2—n 



< cr 2 - n e I( - r) uj{r) for < r < 1. 
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To estimate h\ we use \t(s)\ < ca(s) together with (|73|) and a similar analysis to the above to 
obtain 

(84) |/n(r)| < cr 2 ~ V (r) max(w(r), a(r)) for < r < 1. 
Define £(r) by 

(85) CM = fOT < r < 1. 

h Q {r) 

Using (f8"Tj) - ([HI]) we can estimate \({r)\, |^C' ( r ) I < c max(u;(r), c(r)). To estimate we write 
h (" = h'{ - /i'o'C - 2/i'oC', where 

(86) ftg(r) = (n - 1) r'VM - r 1 " V (r) /'(r), 
and 

^'(r) = r~ V (r) [(n - l)r(r) - rI'(r)T(r) - rr'(r)] 

(87) /- r 

+ r--"£J + (r)(n- 1 + fl(r)) / p n - 1 R(p)B[D 2 v](p)dp - B[D 2 v](r). 

Jo 

The terms /ig'C an d 2/iq£' may be estimated pointwise as before, but h'{(r) involves the term 
B[D 2 v\{r), which cannot be estimated pointwise. However, from (|62p and v (£ X we conclude 
M p (r 2 C",f) < cmax(w(r),ff(r)). Putting this together with the lower order derivatives, we 
obtain the desired estimate (|38p . Summarizing so far, we have found a solution Z of ([35]) in the 
desired form (|36|) . 

Next we need to verify that any strong solution u £ W e ^ p (Bi\{0}) of £u = that satisfies the 
growth estimate (|40]) must be of the form (j4Tj) . To do this, we shall invoke well-known results for 
weighted Sobolev spaces. To begin with, let us introduce the weighted L p -norm on B Q = Bi\{0}: 



(88) HW)=/ \x\f*\u(x)\>dx. 

A)- / : 



To relate this to the M p -norm, notice that 

r. OO r. 

/ \xfP\u(x)\Pdx = Y, \xf p \u(x)\ p dx. 

J0<\x\<l k=1 J2- k <\x\<2- k + 1 

Moreover, there exist constants c\, C2 (independent of k) such that 

Ci2 -k(/3p+n) M ^ 2 -ky < / \x\ 0p \u{x)\ p dx < C2 2- k ^ p+n '> M p (u, 2- k ) p . 

J2- k <\x\<2- k + 1 

Consequently, M p (u, r) < cr a for < r < 1 implies u £ L^(B ) if a + > —n/p, and conversely, 
u £ L p (B a ) implies M p (u, 2~ k ) < c 2^+ fe "/p which implies M p (u, r) < c Q r a for < r < 1 if we 
choose a = —f3 — n/p. We obtain analogous relationships between the M p -norm and the L^-norm 
of the terms \x\^\d a u(x)\ for \a\ < 2. 

Now let us introduce a weighted L p -norm for functions on R™ = R™\{0} with separate weights 
at the origin and infinity: 



(89) 



\u(x)\ p \xf p dx + / \u(x)\ p \xp p dx, 

0<|a:|<l J\x\>l 
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where B c = R"\i?i. We then define the weighted Sobolev space Wp*(M%) to be those functions 
in W^(R") for which the norm 

(90) h\\ w; .» = HM H ^«IIl^ (R? ) 

|a|<2 

is finite. Many authors have used similar weighted Sobolev spaces to study operators like the 
Laplacian on R n , R™, and other noncompact manifolds with conical or cylindrical ends. Using 
the analysis in [TB], [TH] or [IB], for example, it is easily verified that the bounded operator 

(91) A : W|*(R?) - L? +2n+2 (K) 

is Fredholm (finite nullity and finite deficiency) for all values of (3 and 7 except for the values 
— 2 + — + k and — — — k where q = p/(p — 1) and k is any nonnegative integer. In fact, (|9ip is 
an isomorphism for —n/p < /3, 7 < —2 + n/q (recall that we are assuming n > 3, so such /3, 7 
exist). Since we are principally interested in the behavior of functions at the origin, we will fix 
70 G (—n/p,— 2 + n/q). Then 

Moreover, for (3+ G (—2 + n/q, —1 + n/q), we find that (f^Tj) is surjective with a one-dimensional 
nullspace spanned by |x| 2 ~™ G W|^ P 7o (R"); for /?_ G (—n/p — 2,— n/p — 1), we find that ([9l"j) is 
injective with codimension equal ton + 1. 

Next we use perturbation theory (cf. [TO], Ch.IV, Sec. 5) to conclude not only that the operator 

(92) C:W; : P 7 (K)^L^ +2i7+2 (K) 

is Fredholm for exactly the same values of (3 and 7 as for (f9Tj) . but the nullity and deficiency of (f9"Tj) 
and (f9"2")l agree, provided <5 is sufficiently small. So, in addition to the fixed 70 € (—n/p, —2+n/q), 
let us now choose /3 G (—n/p, —2 + n/g), as well as /3_ e (—1 — n/p,— n/p) and /3+ G (—2 + 
n/q, — 1 + n/q). Then we assume that S is so small that 

(93) C : W** m (R; 1 ) - ^ 0+2)70+2 (R™) ^ an isomorphism, 

(94) £: W^ 70 (R?)^^_ +2i7o+2 (R?) is injective with codim = n + 1, 
and 

(95) C : 70 (R o ") - L£ ++2!7o+2 (R^) is surjective with nullity = 1. 

We claim that Z G W^ v (R™). In fact, this is quite simple to check given the explicit formulas 
(|8ip and (|82)l . and the fact that seX, where X has the norm (|53|) . Thus Z is a basis vector for 
the one-dimensional nullspace of (f9"S"]) . 

Now suppose u G W^^(i3i\{0}) satisfies Cu — and the growth estimate (|40|) with e S 
(0,1). Introduce a cut-off function x € Cq°(Bi) equal to 1 on -Bi/ 2 . If we now specify that 
13+ G (-1 - £ + n/q,-l + n/q), then X u G W|; P 7o (R?). Let / = C( X u). Since / = for 
|x| < 1/2 and for \x\ > 1, / G L^ q+2 7o+2 (R"). But (|93|) is an isomorphism, so we can find 
v = Cr x j G W|^ 7o (R"). Now X u - v G VF|; P 7o (R") satisfies £( X u - v) = 0. Since the nullspace 
of ()95p is spanned by Z, there exists a constant C such that x u ~ v — CZ. But this means in 
particular that u = CZ + v for < \x\ < 1/2. 

Now let us describe v asymptotically. Let Co> Ci> • • • ? Cn denote a basis for the cokernel of (IM)) . 
i.e. the C-t are linear functionals on +2 7q+2 (R") that are linearly independent and vanish on 
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the image of ([94]) . Now we want to choose Co, Ci, ... ,C n so that C(x(v — Co — Y^j=i Cj x j)) ^ s 
in the range of (|94|) . i.e. 

n 

(96) (i[£(x(Co + Y,C j x j )]=<; i [£(xv)] fori = 0,...,n. 

3=1 

To be able to solve the linear system for Co, ...,C n , we need to verify that the finite- 
dimensional linear map 

(n n 
Co[£(x(c + c i x i)]> ■ • • > C»[£(x(c + J] Cto) 
3=1 3=1 

is nonsingular. Suppose (Co, . . . , C„) is in the nullspace of (|97| . Then £(x(Cq + Y^=i CjXj] is 
in the range of (|94|) and there exists -0 G (R™) such that 

n 

£(V>)=£(x(C + £C^)), 

3=1 

i.e. £(0 - x(C + E" =1 Cjafj-)) = °- But V> - x(C) + £* =1 € W|^ o (M?) and HMD is an 

isomorphism, so ip — x(Co + Ej=i CjXj). However, \ and x^j are noiin W^ P 7o (]R™), so the only 
way that we can have ip = x(Co + Y^j=i ^j x j) IS t° have Co = C\ — ■ ■ ■ — C n — 0. Consequently, 
the linear map (|97p is nonsingular. 

Thus we can find Co, . . . , C„ and w € WVf (IRq) so that C(w) — C(x(v—Cq—CiXi ■ ■ ■ C n x n )). 
But (|93|) is an isomorphism, so we obtain v = Co + C\X\ + • • ■ + C n x n + w for < \x\ < 1/2. 
This yields (gT|) and we only need to verify M 2:P (w, r) < cr 2 ~ El . But recall that w € W^5, 7o ( K o) 
implies M2 )P (w,r) < cr~P-~ n l v for < r < 1/2, and we can then let /3_ = —2 — n/p + e± for 
any ei G (0, 1) to obtain the desired estimate. □ 

Now let us formulate the result for a general point y € M™ where we do not assume a^- (y) = #y . 
With y fixed and the same conditions ([5]) and on w, we now assume 

(98) sup ||A X - AJ <w(r) for < r < 1, 

|a— j(|=r 

and we want to construct a solution of 

(99) C(x,d x )Z y (x) =0 forxeS £ (y)\{y}, 

for e sufficiently small. (Unlike Theorem 1, in the Corollary below we need to assume that the 
coefficients are real- valued so that we can choose coordinates in which a.y(0) = 6ij.) 

Corollary 2. For n > 3 and p £ (1, oo), fix y £ M 71 and suppose that the constant coefficient 
operator C(y, d x ) is elliptic and the coefficients Uij(x) are bounded, measurable functions satisfying 
i98\) . For e > sufficiently small, there exists a solution of H99\) in the form 

(100) Z y (x)=h y (\A- 1 / 2 (x-y)\)+v(x) 
where h y is defined by 

(101) h y (r) = £ s 1 -" e'*« ds (1 + C y (r)), 
with I y (r) given by \1S\) and 

(102) M 2 , p (( y ,r;y) < c max (uj (r),a(r)) , 
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and v in il00\) satisfies 

(103) M 2 , p (v,r;y) < cr 2 - n e 1 ^ cj(r). 



Moreover, for any u £ W^(B e (y)\{y}) that is a strong solution of C(x,d x )u = in B e (y)\{y} 
subject to the growth condition 

M 2 Ju,r;y)<cr 1 - n+eo 
where e > 0; there exist constants C, Co, C\, . . . , C„ (depending on u) such that 

n 

(104) u[x) (xj ~Vj) + w(x) for < \x - y\ < e, 

3 = 1 

where w satisfies 

(105) M 2tP (w, r; y) < cr 2 ~ £l for any e x > 0. 

If we use integration by parts, we can write the solution in Corollary 1 as 

(106) Z y (x) = (Vfr-^fr-v)) 2 * 1 e /,(V<A,-(x-»),(,-»))) + 

[n 2) 

where Mi t00 (£ y , r;y) < c max(w(r), a(r)) for < r < e. Notice that, if I y (r) has a finite limit as 
r — ► 0, then the leading term in l|106|) is just a constant times F y (x — y), the fundamental solution 
fatC(y,d m ) at y (cf. (0)). 

Proof of Corollary 2. For now we continue to assume y = 0, i.e. 

(107) sup ||A X -A || < w(r), 

\x\—r 

— 1/2 

and let us assume Ao is positive definite. Let B = = A so that BAqB = I. Introduce 
new independent variables x = Bx and the matrix A^ = (ay(x)) = BA^B, which satisfies 
Sij(0) = Sij. Using d/dxi — b^djdx^ = bkid/dxk, we can then write 

(108) C(x,d x ) = ^{x)^- = S«(*)^^ = 

and apply Theorem 1 to £ in the coordinates x. We conclude that for e > sufficiently small, 
there is a solution Z of £(£, d x )Z(x) = for < |x| < e, of the form Z(x) — h(\x\) + v(x) where 
M 2 , p (v,r) < cr 2 ~ n e I( - r ^>u(r) and h(r) is of the form (J37J) with 

(109) /(r) = ^_/ ftrA 2 -n(A^,2)|2|- 2W/: 



15" 



r<|z|<e 



2 



Expressed in terms of the original variables x, we obtain .Z(x) = Z(x) = h(\Hx\) + v(Bx) that 
satisfies C(x, d x )Z(x) = for < |Bx| < e; but choosing eo sufficiently small, we conclude that 
C(x, d x )Z{x) = for < \x\ < e . 

y i 2 ~ 

Finally, if y is a general point in R", then let x = B(x — y) with B = Aj, ' and let A x = 
(dij(x)) = BA^B; since x — corresponds to x = y, we have a,ij(0) = Sij and we can apply 
Theorem 1 to C(x,d x ) = C(x,d x ). We obtain the solution Z(x) = h(\B(x — y)\) + u(B(x — y)) 
where h(r) involves I(r) as in (|109[) . To transform p09[) to the original variables, replace Aj by 
A z and every other occurrence of z by z — y; we find that (|109[) is of the desired form (fT3")) , so 
we may let h y (r) = h(r). Moreover, since v satisfies M 2 ^ p (v,r) < cr 2 ~"e / ^ r ^(r), it is clear that 
v(x) = v(B(x — y)) satisfies M 2tP (v, r; y) < cr 2 -" e '('') w ( r ) j completing the proof. □ 
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4. Finding the Constant C y so that CZ y {x) = C y 5(x — y) in B e (y) 

Let us now begin to discuss the role that Z y (x) plays in finding the fundamental solution. As 
before, we first consider y = with eijj (0) = 5^ and then use a change of variables to consider a 
general fixed y £ M n ; as in the proof of Theorem 1 , we shall assume e = 1 . We first want to see 
whether the function Z(x) found in Theorem 1 satisfies 

(110) - C{x, d x )Z{x) = C S(x) forxGBi(O) 
for some constant Co- 

It is not immediately clear how the left-hand side of (|110| is defined. Recall from the proof of 
Theorem 1 the decomposition Z(x) = h(\x\) + v(x), where h is given by ([80| and M2 yP (v,r) < 
cr 2 ~ n e 1 ^ ui(r). We can easily calculate didjZ and show that for any fi > there is a constant 
Cfj, so that 

(111) M p {d i d j Z,r)<G ll r- n -^. 

In fact, since the are bounded functions, we conclude that for each i,j the function 

(112) F lj (x) = a ij (x)d i d j Z(x) for x^O 
satisfies 

(113) M p {F tjl r) < C p r- n ->* for < r < 1. 

But this estimate implies that Fij can be regularized at x = to give a distribution Tij, in 
particular as a continuous linear functional on the space of A-H61der continuous functions of 
compact support in U = Bi(0): 



(114) (Fi3,<t>)= F tJ (x)[0(x)-mx(\x\)]dx 

J\x\<l 

where x( r ) i s a smooth cut-off function which is identically 1 near r = but vanishes for r > 1/2. 
(Since \4>{x) — 0(O)x(|a;|)| < C*|a;| A , by choosing < fi < A we see that the integral in (|114[) 
converges.) Now let us take the sum over all 

(115) r = Y,^ir 

Then is a regularization of C{x, d x )Z(x) — 0, and so it vanishes: 

(Jb, 4>) = Yjfa,*) = / E F «( x ) ~ </>(°)x(\x\)] dx = o. 

Of course, regularization effects the distribution only at x = 0, so if we can interpret the expression 

(116) T = C{x,d x )Z{x), 

as a distribution, then T is supported only at x — 0. As such, it is a linear combination of the 
delta distribution and its derivatives. But since J- is a continuous linear functional on A-H61der 
functions with A € (0, 1), it must only involve the delta distribution itself, i.e. (II 10[) must hold for 
some constant Co- But we still have two problems: (i) how is T itself defined as a distribution? 
and (ii) how do we calculate the constant Co? 

The difficulty in defining T as a distribution in U is a consequence of the lack of regularity of 
the aij. In particular, there is no difficulty in defining the distributional derivatives of Z: 

(117) (d l d J Z,<j ) ) = - [ djZ(x)di<p(x)dx 

Ju 
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for (f> 6 Cq(U), since the integral on the right converges. To handle the ay, we can write 
(118) (CZ, <(,)=[ [(ay (a;) - 5 ij )d i d j Z(x) (f>{x) - d i Z(x)d i ^(x)] dx. 



Of course, the integral in (|118|> is actually an improper integral due to the singularity of didjZ 
at x = 0; but provided this integral converges, we conclude that — CZ = CqS in U, and we can 
calculate Co by 



(119) 



Cn = lim 

£->0 



|x|<l 



[(-aij(x) + Sij)didjZ(x) (f> s (\x\) + diZ(x)di4> £ (\x\)}dx, 



where <^> e (|:r|) = with the cut-off function \ introduced above; for these purposes, we are 

able to assume <p(x) = <fi(\x\) since we can write 4>(x) — 4>o(\x\)-\-if)i(x) with |0i(a;)| + |a;| |V0i(x)| < 
c\x\ for \x\ < 1, which shows that (CZ,<px) is well-defined and cj>i contributes nothing to Co. We 
shall now prove the following. 

Theorem 2. Suppose the conditions of Theorem 1 hold and Z is the junction found there. If 
1(0) = lim r _^o I( r ) exists and is finite then the integral in M18\) converges and we can calculate 
Co = \S n ~ 1 \e 1 ^ . If I(r) — > — oo as r — > 0, then the integral in U18\) converges and we can 
calculate Co = 0. 

Proof of Theorem 2. Recall from the proof of Theorem 1 the decomposition Z(x) = /i(|a;|) + 
v(x), where h is given by and v satisfies Since we always assume that I(r) is bounded 

above, we obtain from ([39")) 

(120) M 2 , p (v,r) < cr 2 ~ n uo(r) for < r < 1. 

We shall separately consider the roles of v and h. In the estimates below, f\ x \ <£ should actually 
be considered as an improper integral lim^o X ) <| a |< e ' Du t, we avoid such cumbersome notation. 
For v, we use (I5TT) and (|120[) to conclude 



diV di6 F dx 



|x|<e 



< ce 



ui(r) dr < cuj(e) — > 0. 



and 



|x|<E 



<5y) « 4>s dx 



< c 



/ u (r)—= co- 
Jo r 



(s) - 0. 



This shows that v makes no contribution to the value of Co- 

To determine the effect of h, let us write h(r) = h 2 (r) + h 3 (r) + c where 

„1— n ) 



(121) 



h 2 (r) =ci / s L ~ n E + (s)ds 



with ci chosen as in the proof of Theorem 1 so that C\ E + (r) 



oil) as r 



0, and 



(122) 



h 3 (r) 



l E + (s) 



l E_(p)B[D 2 v}(p)dp. 



(Notice that h 2 and /13 differ slightly from ((81]) and (|82[) .l Let us consider /13 first: 

9^3 = -x ir - n E + (r) r p n - 1 E_(p)B[D 2 v}( P )dp 
Jo 



and 



d l d j h 3 = 



l E_ 



R(r)) / p»- 1 S_(p)S[£> 2 «](p)dp 



B[L> 2 u](r). 
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The calculation of dihz combined with (f73|) and the assumption that I(r) is bounded shows 



(123) 



dih 



i ri3 <->i<Pe 



dx 



\x\<e 



< ce^ 1 [ a(r)dr < ccr(e) -> 0. 
Jo 



On the other hand, the "worst" term in (a,j — 6ij)didjh,3 is 



W(x) = -r- n E + (r) (a tl - n ^fl ( 



Qij XiXj 



r * \ r 

Using the fact that cf> e only depends on r and observing that 



an - n- 



l)R) I p^E^BlD^dp. 
l)R)dB = IS^lRir), 



we can calculate 



R(r 



(W,cj )£ ) = -\S n - 1 \ E+(r)-^- Mr) p n - l E.(p)B[D 2 v](p)dpdr 



l^"- 1 ! / (E+)'(r)Mr) / P^E^BiD^dpdr 



Now we can integrate by parts to obtain 



{W,4>e) = \S n - 1 \ / E + (r)<j>' e (r) p^E^B^v^dpdr 



+ J ^ e (r)r B_1 B[D a t;](r)(ir 
But, using \<f/ e \ < c£ _1 and \E + \ < c, we find 

E + (r) ft (r) f p"- 1 ^ (p)£?[£> 2 v] (p) dp dr 



< C£ 



and 



e (r) r n_1 -B[£) 2 u](r)dr 



<r(r) cir < ca(e) — ► 0, 
< cer(e) -> 0. 



We conclude that /13 makes no contribution to Co- 

Now consider /i2- We compute dih 2 = —c\Xir~ n E + (r) and 



so 



Wa = - Cl r-» £7+(r) (% - - ^fl(r)) , 

(a,j - 5 ij )did j h 2 = -ci r~™ £+(r) (a« - b "'^'^ - 



Notice that 



M<6 



r""B + (r) (o«-n^^-( 



1)R 



-1 R X (-)dx 



|.V" / £ + ( r )^ x (-)dr, 
r e 



so we now may compute 



-.r/,,.0=,,|,s"'-'| / (S + (r)^x(^)-^ + (r)[x(;) 



dr 



-ci|^ 



S+(r)x(- 



dr. 



Now, if 1(0) exists and is finite, then CiE + (0) = e I(0 \ and we conclude -{Ch 2 , <f> e ) = \S n - 1 \e I ^°\ 
On the other hand, if 7(r) — * —00 as r — * 0, then E + (r) — > 0, so C = 0. This completes the 
proof of Theorem 2. □ 
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Now let us consider a general fixed y £ W 1 and try to find the constant C v so that 
(124) - C(x, d x )Z y {x) = C y 5(x) for x E B x {y). 

We replace (|118|) by 



(125) (£Z V ,4>)= / [{aij{x) - a lj {y))d l d j Z y {x)(j){x) - a ij {y)diZ y (x)d j 4>{x)] dx. 

Ju 

Corollary 3. Suppose the conditions of Corollary 1 hold and Z y is the function found there. 
If I y (0) = linv^o I y (r) exists and is finite then the integral in \125\) converges and Z y satisfies 
{12$ with C y given by 

(126) C y = l^lv/detA^W. 

// I y (r) — > — oo as r — > then the integral in H125\) converges and Z y satisfies {12$ with C y = 0. 

Proof of Corollary 3. We need only show {—CZ V ,4>) = |S , ™ _1 |(detA 2/ ) 1 / 2 e 7 ''^ ^(/)(y) for some 
<j> G C^° {B £y (y)). Let us recall the change of coordinates used in the proof of Theorem 1, namely 

— 1/2 ~ 

x = B(x — y) where B = A y , and let 4>(x) = 4>(x); then 

C(x, d x )Z y (x)4>(x) dx = -(detA^) 1 / 2 / £(x, dx)Z (x)<j>(x) dx. 



But Theorem 2 implies -{£.Z Q ,<j>) = IS" 1-1 ^ 7 "^^). Since ^(0) = 4>(y), we obtain the desired 
result. □ 

5. Constructing the Fundamental Solution 

Now we are in a position to construct the fundamental solution in a bounded open set U C K n 
when dij € C U (U), assuming that lo satisfies the square-Dini condition ([9]) and for every y G U 
we know that I y (0) exists and is finite. Given the results in the preceding section, it is natural to 
seek the fundamental solution in the form F(x,y) = Z y (x)/C y +v(x, y) where C(x, d x )v(x, y) = 0. 
In fact, it will be convenient to construct F(x,y) as the Green's function G(x,y), in the sense of 
||5J|, for a smooth, bounded domain V that contains U. 

But first let us observe that the additional assumption that I y (0) exists and is finite allows us 
to improve the asymptotic description of Z y (x). In fact, let us fix y = with 0^(0) = <5.y, and 
assume 

(127) |/(r)- I(0)| <0(r) 

where 0{r) is a positive, nondecreasing function for < r < 1 such that 6(0) = 0; as with u> we 
additionally assume that for some v € (0, 1) we have 

(128) 9(r) r~ 1+lJ is nonincreasing for < r < 1. 
(If u(r) satisfies the Dini condition, then we can take 8(r) — u>(r).) 

Lemma 1. Under the conditions of Theorem 1, let us additionally assume that {127\ l and A128\) 
hold. Then the solution Z found in Theorem 1 satisfies 

I |2-n 1(0) 

(129) Z(x) = U (1 + ^)), 

n — 2 

where for any p £ (1, oo) we can estimate £ by 

(130) Af 2 , p (£, r) <c max(w(r), cr(r), 0(r)) as r -> 0, 
where c depends only on lo, 9, n, and 1(0). 
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Proof of Lemma 1. In the decomposition Z(x) = h(\x\) + v(x) as in (|36|) . it is elementary to 
use (|39j) with (|127|) to show that (|130l) applies to £ = c |x|" -2 v, so we focus on h. Recall the 
decomposition h = ho + hi using (151]) and (|82p. Let us recall /iq from and use (|127[) to 
estimate e / ^ r - ) < (1 + 2#(r))e 7 (°) for r sufficiently small; also recalling |/'(r)| < 2(n — l)u)(r)/r, 
we obtain 

\h' \r) - (n- l)r-"e 7(0) | < C r" max(w(r), 0(r)); 

here (and henceforth) c depends only on w, 9, n, and 1(0). Using (|34[) and (|128p , we can integrate 
this twice to obtain 

r 2-n e i(0) 

< cr 2_n max(w(r), 0(r)). 



h (r) 



n-2 



In fact, this argument has shown that £o( r ) = i n ~ 2) r ™ 2 e /< "°' ) ^o( r ) — 1 satisfies the pointwise 
estimate 

M 2 ,oc{S,o,r) < c max(w(r),0(r)). 

Similarly, let us recall h'[ from (|87[k we can estimate most of the terms pointwise, but the term 
B[D 2 v\(r) can only be estimated in M p , so £i = (n — 2) r™ -2 e _/ ^ ^ ft-i(r) satisfies 

M 2 , p (£i,r) < c max(w(r),cr(r)). 

These may be combined to yield (|130[) . so the Lemma is proved. □ 

Now we are ready to construct our fundamental solution in a bounded open set. 

Theorem 3. Suppose (QJ) is uniformly elliptic in a bounded open set U C M™ for n > 3, where 
the coefficients aij are continuous functions with uniform modulus of continuity to satisfying 
// I y (0) exists at every y £ U and 

(131) \I y (r) - I v (0)\ < 0(r) for all y E U, 

where 9 is a positive, nondecreasing function for < r < 1 with 9(0) = and il28\) , then there is a 
function F(x,y) satisfying for x,y € U; moreover, F(x,y) admits the asymptotic description 

(A~ 1 (x — w), (x ~ y))^~ 

< i32 » *•"■■»= il+H(xM - 

where for any p £ (1, oo) and any compact set K C U we have 

(133) M 2 , P (H(;y),r;y)<cmax(w(r),o-(r),6(r)) as r -> 0, 
wii/i constant c independent of y £ K . 

Proof of Theorem 3. For each y £ U, denote the e in Corollary 1 by e y . Note that the size of 
e y depends on the behavior of the coefficients a,y near y through their ellipticity (i.e. the norm of 

— 1/2 

A y ) and their continuity (i.e. u>). Since the ellipticity and modulus of continuity are uniform 
on U, we can find e > that is independent of y £ U. In fact, if we choose a smooth, bounded 
domain V with V D U and dist(J7, dV) > e, then we can extend the coefficients <Zy to V in such 
as way as to maintain the uniform ellipticity as well as the modulus of continuity u>(r) (cf. [20] ). 
Thus for all y £ U we can construct Z y (x) in B £ (y)\{y} . In fact, repeating this argument with 
a neighborhood of V, we may assume that Z y (x) is defined for all y £ V. 

For each y £ V let us use Corollary 2 to calculate C y > 0, and conclude that — C(x, d x )Z y (x) /C y - 
S(x — y) for all x,y £ V with \x — y\ < e. We shall construct G(x,y) as the Green's function for 
C in V. For fixed y £ V, let us introduce a smooth cut-off function r/ y (r) satisfying rj y (r) = 1 for 
sufficiently small r > but rj v (\ ■ —y\) has compact support in V. Then let us define 

(134) G(x,y) = riy(\x-y\)Z y (x)/Cy +v(x,y), 
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where v{x,y) is to be determined. But if we apply —C(x,d x ) to G(x,y) we obtain 

-C(x, d x )G(x, y) = 6{x - y) + ip(x, y) - C(x, d x )v(x, y), 
where ip(-,y) £ L P (V). So, for fixed y £ V, consider the Dirichlet problem for v: 

C(x,d x )v(x,y) =ip(x,y) for x £ V, 
( ' v(x,y)=0 forxedV. 

It is well-known (cf. Theorem 9.15 in [5]) that (|135p has a unique solution v(-,y) S W 2 ' P (V) fl 
Wp' p (y), so using this v(x 1 y) in (|134[1 not only ensures that ([8]) holds for x,y £ V (and hence 
for x,y £ U), but the following: for any / £ C(V), u(x) = — J v G(x,y)f(y) dy satisfies Cu = f 
in V and u = on dV. Now if we pick <f> £ C 2 (V) with <fr = on dV and let / = C(f>, then 
uniqueness of the solution shows u = <p in V, i.e. |8]) holds. In other words, the G(x, y) that we 
have constructed is just the Green's function for C in V; in particular, ([2]) holds. 
Now from (|134|) and Corollaries 1 and 2 we see that 

(136) G(x, y) = {A v 1{X ~ y) ' {X ~ e ^(v / (A y - 1 (x-,),( 3: -,))) (1 + Ux)) 

where Mi j00 (^j / , r; y) < c max(w(r), <r(r)) for < r < e and J a (r) = I y (r) — I y (0) can also be 
written as 

(137) -I*-/" LaA- 1 ) n (A z A- 1/2 (z~y),A^ 1/2 (z-y)) \ 

However, if we choose coordinates in which y — and ay (0) = dij, then we may apply Lemma 1 
to absorb the exponential term in (|136[) into 1 + £ y and obtain 

(A~ 1 (a; — y). (x — y))^~ 

where M 2iP (^ y ,r;y) < c a max(w(r), c r (r), 6*(r)) for < 7- < e with c y depending on w, 0, n, and 
1^(0). But if we select a compact subset K C U, then c y may be taken independent of y £ K, so 
we can replace £, y (x) by H(x, y) and obtain (|132p . (|133[) as |a; — y| — » 0. 

Letting F(x, y) = G(x, y) for x,y £ U C V , we have our fundamental solution in U. □ 
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